Quantum Bounds on Bell inequalities 
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We have determined the maximum quantum violation of 241 tight bipartite Bell inequalities with 
up to five two-outcome measurement settings per party by constructing the appropriate measurement 
operators in up to six-dimensional complex and eight-dimensional real component Hilbert spaces 
using numerical optimization. Out of these inequalities 129 has been introduced here. In 43 cases 
higher dimensional component spaces gave larger violation than qubits, and in 3 occasions the 
maximum was achieved with six-dimensional spaces. We have also calculated upper bounds on 
these Bell inequalities using a method proposed recently. For all but 20 inequalities the best solution 
found matched the upper bound. Surprisingly, the simplest inequality of the set examined, with 
only three measurement settings per party, was not among them, despite the high dimensionality of 
the Hilbert space considered. We also computed detection threshold efficiencies for the maximally 
entangled qubit pair. These could be lowered in several instances if degenerate measurements were 
also allowed. 

PACS numbers: 03.65.Ud, 03.67.-a 



I. INTRODUCTION 

Suppose two classical systems which are separated 
from each other. Let us make some local measurements 
on them. Then the results of these measurements may 
be correlated, which may be explained by shared ran- 
domness experienced in the past. For a given number of 
possible measurement choices (inputs) and results (out- 
puts) the set of correlations forms a polytope whose finite 
number of vertices correspond to all the deterministic as- 
signments of outputs to inputs. The facets of the poly- 
tope, which form the boundary of the classical region, 
correspond to tight Bell inequalities P, Thus Bell 
inequalities has no a priori relation to quantum physics. 
However, quantum physics violates them, which has been 
verified experimentally in numerous occasions up to some 
technical loopholes (see e.g., Q). Indeed, this makes Bell 
inequalities very interesting. 

On the other hand, one may ask what is the achievable 
set of correlations if one allows the two parties to share 
quantum resources as well over shared randomness. This 
is a convex set, such as for the classical case, however it 
cannot be described by a finite number of extreme points 
Nevertheless, one can construct so called quantum- 
Bell inequalities which bound the correlations achievable 
by quantum physics (see also item 5 in Sec. III.C 

of Ref. 0). A simple way to form quantum-Bell inequal- 
ities is to use the coefficients of known Bell inequalities, 
and determine the maximum value one can get by per- 
forming measurements on quantum systems of arbitrary 
dimensions. 

The simplest result in this respect is the Tsirelson 
bound Q stating that quantum correlations cannot vio- 
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late the CHSH inequality by more than (^/2 - l)/2. 
Note, however that the requirement of non-signaling 
alone allows a higher bound of 1 (the maximum 

allowed value in a local classical theory is 0). In the 
present paper we stick to the bipartite two-outcome sce- 
nario and we consider inequalities with more than two 
measurement settings having nontrivial local marginals 
as well. 

In particular, in this paper we test 241 bipartite tight 
Bell inequalities with two-valued observablcs with up to 
five settings. These Bell inequalities have been collected 
from various sources: from the list of Ref. [lH contain- 
ing 89 inequalities labelled by Ai, i = 3,..., 89, from 
the list of Ref. comprising 31 inequalities with up to 
four settings per party. 18 of them, labelled by 14422 , * ^ 
3, . . . , 20 are introduced by that paper. The simplest non- 
trivial inequality 13322 beside the CHSH inequality was 
presented first by Froissart Il3l i n 1981 and more recently 
it was reintroduced by Sliwa |14l | and by Collins and Gisin 
[Tsl . The rest of them, namely, ll!^2 ^^'^ -^4422 were pre- 
sented in Ref. 14422 0> ^5,6 and AIIi^2 in [13 j 
while AS I ^2 in Q- These lists were also considered by us 
previously [l^. Finally, a list of 129 tight bipartite Bell 
inequalities with four settings per party is constructed 
in the present paper (labelled by J4422, « = 1, ■ • ■ j 129). 
We test numerically these Bell inequalities and give the 
bound on them achievable by quantum systems, thereby 
generating quantum-Bell inequalities. On one hand, we 
used numerical optimization in higher dimensional (up to 
six complex and eight real) component Hilbert spaces to 
calculate a lower bound on maximal Bell violation. On 
the other hand, the hierarchy of tests of Refs. [l^, [l^l 
allows us to give upper bounds on the violation achiev- 
able by quantum systems. We found in a surprisingly 
many cases that the lower and upper bounds coincide 
indicating that the quantum bound we have found is 
tight. Numerics suggests that in some cases even five- 
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dimensional local Hilbcrt spaces (i.e., 25-dimensionaI full 
Hilbert space) were not enough to achieve the maximum 
quantum value. This result can be interpreted through 
the recent concept of dimension witness 2l[ , which is a 
certain kind of generalization of Bell inequalities which 
enables to measure the dimension of the Hilbert space 
(see also Refs. [l^, [2^ for recent works using this con- 
cept). In particular, a d-dimensional Hilbcrt space wit- 
ness allows to distinguish the strength of correlations of 
measurement outcomes on members of composite sys- 
tems that can be obtained in c?-dimensional component 
Hilbert spaces from the ones which are achievable only 
in higher dimensions. Our result thus indicates a five- 
dimensional witness, which improves on our recent result 
p^ . where we found three-dimensional witness. We also 
made a study, given these tight Bell inequalities, from a 
different perspective. Namely, we computed the thresh- 
old detection efficiency for a maximally entangled qubit 
pair required to close the detection loophole in the sym- 
metric (Alice and Bob having the same detection effi- 
ciency) and in the asymmetric cases (Alice has perfect 
detection efficiency) by considering degenerate measure- 
ments (having deterministic outcomes) as well. Surpris- 
ingly, in almost all cases degenerate measurements were 
required to achieve the lowest threshold efficiency in both 
the symmetric and asymmetric situations. However, let 
us note, that these cases with degenerate measurements 
could be traced back to ones with non-degenerate mea- 
surements, which achieve the same threshold value. In 
the asymmetric case we could improve on the best known 
results from the literature [l^, Hj] , by carrying out only 
non-degenerate measurements. 

The rest of the paper is organized as follows. In SecHIl 
we describe methods to find tight two-outcome bipar- 
tite Bell inequalities, and we show the 129 inequivalent 
examples we have got this way. For all of them the clas- 
sical limit is zero, and we have given a symmetric form, 
whenever there exists one. In Sec. Illll it is explained how 
we constructed the maximum quantum violating solu- 
tion in real and complex six-, and real eight-dimensional 
component Hilbcrt spaces. In Sec . IIV I we present briefly 
the method of Navascues et al. 

mm 

to obtain up- 
per limit for the maximum quantum violation of two- 
party two-outcome Bell inequalities. In Sec. |V]we give 
the maximum quantum value we have derived, the size of 
the Hilbert spaces required, and the dimensionalities of 
the projectors used as the measurement operators for all 
the 241 Bell inequalities considered. We also show when 
this value reaches the upper bound. For those cases the 
value given is the absolute quantum bound. We discuss 
the results obtained. In Section IVII the minimal detec- 
tion efhciencies are computed for a maximally entangled 
pair of qubits considering degenerate measurements as 
well. More details of the results including the matrices 
of the optimal operators and the Schmidt coefficients of 
the states are presented in the web site [s^ ■ Section IVIII 
summarizes the results of the paper. 



II. METHODS FOR GENERATING TIGHT 
BELL INEQUALITIES 

We consider a typical experiment to test for correla- 
tions. Imagine a system, cither quantum or classical, 
which is composed of two subsystems, A and B. Suppose 
that Alice has choices of two-valued measurements 
on A, while Bob has tub choices of two-valued measure- 
ments on B. Let us refer to this situation with the no- 
tation myiTOB22. The outcome of the measurements is 
one of two values labelled by the symbols 1 and 0. The 
experiment is repeated many times with different choices 
of measurement settings so as to get accurate probabili- 
ties. The result of such a correlation experiment can be 
thought of as a vector p in a dimensional space, where 
d = m.A + TUB + ruATriB- The components PAi^ Pb, and 
PAiBj (for 1 < i < ruA and 1 < j < ras) of p represent 
the probability that Alice gets 1 for measurement setting 
i, Bob gets 1 for measurement setting j, and that both 
Alice and Bob get 1 for setting i and j on the respective 
sides. 

In a quantum correlation experiment Alice and Bob 
share a quantum state p, and the probabilities can be 
expressed by the formulae 



PA, = Tr[p(A, ® /b)] , 
PB, = Jr[p{lB<E)B/)] , 
PA^B, = Tr[piA,(g) Bj)] , 



(1) 



where I a (Ib) denotes the unity operator acting on the 
Hilbert space of Alice (Bob), while Ai and Bj are the 
{1,0} valued observables (i.e., projectors) measured by 
Alice and Bob, respectively. Note, that for our case of 
binary outcomes it has been shown [25j that projective 
measurements can reproduce the set of all quantum cor- 
relations. 

On the other hand, the result of a classical correlation 
experiment must correspond to a probability distribu- 
tion represented as a convex sum over all deterministic 
configurations. A deterministic configuration assigns the 
outcomes {1,0} to each of the tua + ttib measurements, 
in which case the joint probabilities factor into two local 
probability distributions, PAiBj =PAiPBj- Thus, the set 
of vectors p which are possible results of a correlation ex- 
periment forms a, d = {uia -I- mp -I- my^ms )-dimensional 
polytopc P [l^. Froissart [ij has shown that in the 
general mA'mB'22 case the set of all n = 2'"-^+'"b points 
Pu, V = l,...,n corresponding to deterministic config- 
urations are actually the extremal vertices of this d- 
dimensional polytopc P. The experimental result has 
a locally causal model if and only if the corresponding 
point is located inside the polytope P. A Bell inequality 
is a linear inequality satisfied by all the points in such 
a polytopc. A tight Bell inequality, which is the most 
useful to detect non-local correlations, defines a facet of 
this polytope. In our particular formulation we represent 
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a Bell inequality in the form 

VIA inB jyiA niB 

^ bA.PA, + ^ bBjPBj + X! X! bAiBjPA.Bj < bo, (2) 
i=l j = l i=l j = l 

where the sums run over all measurements of Alice and 
Bob, and the coefficients 5's are suitably chosen real num- 
bers. The simplest nontrivial tight Bell inequality (with 
the smallest values of uia = 2 and ms = 2) is the CHSH 
inequality 0] 

-PAi -PBi +PA1B1 +PA1B2 +PA2B1 -PA2B2 < 0. (3) 

Also, note that for a given Bell inequality bp < bo, with 
p G and for a d-dimensional poly tope P, the face 
represented by the inequality is a facet of that polytope 
P if and only if the rank of the matrix containing the 
vertices py of polytope P which saturate the inequality 
is d — 1. This fact enables us to test tightness of Bell 
inequalities. 

Given the vertices py of polytope P in case mAmB22, 
our task is to find some facets, which defines tight Bell 
inequalities. A deterministic algorithm (see for instance 
the package cdd [131 ), which explores systematically all 
the facets, runs too slowly already for the 4422 case to 
complete this task. This fact was also mentioned by 
Collins and Gisin in Ref. [l^. In fact, it has been shown 
by Pitowsky that this problem for the general mAmB22 
case is NP-complete [20|. Therefore, we turn to heuris- 
tic algorithms and provide two methods which enable us 
to produce tight two-outcome bipartite Bell inequalities 
with mA = mB = 4: measurement settings. The first 
one uses some portions from Seidel's shelling algorithm 
[sot , while the second one is based on a reduction of the 
polytope P to smaller ones. Actually, the second method 
proves to be useful for a larger number of measurement 
settings as well. For instance, we found that it also works 
for the 5522 case. 



A. Algorithm: shelling of a Bell polytope 

In this subsection we suitably modify a part of Sei- 
del's shelhng algorithm [s^]. Consider a d-dimensional 
polytope P corresponding to the case mA'mB'2.2 (i.e., a 
Bell polytope), and let S denote the set of vectors Pv 
corresponding to deterministic configurations as defined 
earlier. Imagine travelling along a directed straight line 
L that intersects the interior of P. Start at a point on 
L that is in the interior of P and move along L. Con- 
tinue moving until you pass through some facet F and 
then leave the interior of P. This facet F is the first in 
the shelling order. Continue moving away from P along 
L and more and more facets will become visible. The 
order in which these facets appear is the shelling order. 
This algorithm is deterministic in essence, which makes 
a shelling of the polytope by travelling along a line and 
determine the facets of the convex hull as they become 



visible. However, this algorithm runs too long to com- 
plete the task in a reasonable time for polytope P of the 
case 4422. Therefore, a heuristic algorithm is given which 
uses a part of Seidel's algorithm: 

Suppose we know a priori a single facet F of the poly- 
tope P defined by all vertices Py £ Sp which lie on the 
supporting hyperplane of this facet F. Define a by adding 
up these vectors p^ . This sets the direction of the shelling 
line L, parameterized by x{t) = —a/t, —00 < t < +00. 
Linear programming can be used to determine for every 
Pv G S the first facet in the shelling that contains py. 
Thus for each py G S one obtains a tight Bell inequal- 
ity defined by a facet of P. Let us collect these facets, 
which serve as new inputs to the algorithm. This rou- 
tine is iterated as long as no more new Bell inequalities 
are found. However, several of them may turn out to be 
equivalent, i.e., symmetric under party exchange, observ- 
able exchange or relabelling of outcomes (see e.g., [3l|). 
Thus one has to take care of selecting the inequivalent 
Bell inequalities. 

B. Algorithm: slicing a Bell polytope 

Similarly as in III Al we start from an a priori known 
facet _F of a d-dimensional polytope P corresponding to 
a tight Bell inequality, bp < bo, where p G W^. In this 
case for all vertices py G S corresponding to deterministic 
configurations the inequality holds by definition. Now let 
us decrease the value of bo to a smaller value bo ■ Thereby 
we obtain two groups of vertices, one (5') for which bp> 
6q holds and one {S") for which it does not hold. 

The convex hull of these vertices form polytopes, des- 
ignated by P' and P" , respectively. By setting an appro- 
priate value for 6q wc may get a polytope P' in the same 
dimension d as for the polytope P. Now we take as an 
input these vertices to the cdd program [23] which gener- 
ate all the facets of this reduced polytope P' (assuming 
that the number of vertices in P' is small enough to be 
handled by the program). However, not all of the gen- 
erated facets will be the facets of the original polytope 
P as well. Thus one must sort out the valid faces which 
fulfill b'py < b'o for all py G S, with b' and 6g defining a 
facet of P'. Such as in the algorithm III Al these facets, 
which define tight Bell inequalities, form the input to the 
new cycle of computations. The computation halts only 
when no more new Bell inequalities can be produced this 
way. 

C. Comments specific to the 4422 Bell inequalities 

A crucial point of both algorithms is the right choice of 
the initial facet. A facet in the d-dimensional polytope 
P is called degenerate if it is created by more than d 
points. The problem lies in the fact that a Bell polytope 
P in general, and also in the particular case of 4422 is 
highly degenerate. However, it is not favorable to start 
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TABLE I: Tight bipartite Bell inequalities with four two-outcome measurement settings per party, generated with the methods 
of Section ini The coefficients, positive or negative single-digit integers, are shown in the order of 6ai, bA2, bA^, bA^^, bsi, &S21 

6S4, &A1S1, &A1S2, ^AiSa, &A1S4, ^AaSi, bA2B2, '^A2S4i ^AjSi, bAsB^, bA^B^, &A3B4, 6A4B1, bA^B^, bA^B^, bA^B^ 

(see Eq. Q)- For each inequality only one of the equivalent forms are shown. The form chosen has classical value zero, and 
whenever possible, it is symmetric. 
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00-2-50-4-6-2-3212-123-42-3423522 

0-2-4-3-5-20032-1-222-2212314-421 

-3-100-3-1003112111-211-2-12-2-1-1 

0-3-4-4-5-20042-1-322-2312314-422 

00-3-60-4-7-3-4213-224-43-4534622 

-3-30-2-3-30-2412-2112322-42-232-1 

00-1-40-2-4-1-2111-112-21-2213321 

00-1-40-2-4-2-3112-112-21-2213322 

-200-4-200-422-322-3-11-3-1-132133 

-8-20-2-8-20-234534-62252-1-432-44 

-10-2-3-10-2-3-21121-1-321-3322221 

0-1-3-1-3-10-121-1-111-2112212-211 

-2-300-3-10-2211-2122211-212-201 

-5-2-10-5-2-10123322-323-301321-3 

-60-3-3-1-4-104222121-332-21-5413 

-50-2-1-2-4-303233222-422-31-3311 

-40-1-1-1-3-202222121-311-21-2211 

-5-1-2-10-5-102322131-313-31-3311 

-3-1-1-10-3001211021-212-20-2201 

-1-4-1-2-6-200321-2122232-314-411 

-2-50-2-6-20-1322-3123332-314-411 

-20-4-20-30012-1102-1-22220-3311 

-30-2-1-2-2-302122212-321-21-2211 

-10-6-5-1-4-1012-21121-34222-5433 

-200-5-30-2-2-32121-1-112-3214411 

-20-2-3-20-2-34-412-4-21311-122321 

-1-10-4-20-3-1-1111112-31-3212222 

00-30-30-2-4-3-1121-1-1233222-310 

00-50-5-2-1-321-1-211-3234432-311 

00-2-50-2-4-1-3121012-32-2213322 



the algorithm III Bl from a highly degenerate facet, since 
in this case the number of vertices of P' which are fed 
in the cdd program can easily become larger than the 
number the program can handle. In our specific case we 
used Brunner and Gisin's l\^22 ^^^^ inequality as an input 
facet to the algorithm and the 129 tight Bell inequalities 
were obtained using both methods HT^ and [TTbI How- 
ever, we mention that we did not perform an exhaustive 
search, since it would have been time consuming. We 
suspect that the list of 26 inequalities of Ref. [I4I plus 



our 129 introduced ones do not form the complete set of 
all 4422- type of tight Bell inequalities. This is implied 
by the fact, that our set of 129 inequalities covered only 
23 from the list of 26 inequalities of Brunner and Gisin. 
The coefficients defining the 129 inequalities are shown 
in Table U 
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III. LOWER BOUNDS 

The quantum value of a Bell expression is the expec- 
tation value of the operator 

rriA niB 

bA, {A,®Ib) + Y^ he, {Ib ® Bj) 

1=1 'j=i 

+ Y.Y.^A'B,{A,®Bj). (4) 

According to Eqs. ([T|) and ^ this is the value one gets 
for the left hand side of Eq. ([2|) with measurements on a 
quantum system. This can be larger then the maximum 
classical value 6o, the quantum violation is the difference 
of the quantum value and 6o- For the maximum quantum 
violation it is enough to consider pure states. 

We have given lower limits for the maximum vi- 
olation of Bell inequalities by explicitly constructing 
measurement operators in Alice's and Bob's compo- 
nent Hilbcrt spaces and the appropriate state us- 
ing numerical optimization. We have made calcula- 
tions with six-dimensional real and complex spaces and 
eight-dimensional real spaces, taking three- and four- 
dimensional projection operators as measurement oper- 
ators, respectively. The results of such a calculation 
may correspond to solutions in component spaces of 
lower dimensionality. The six-dimensional case incor- 
porates the five-dimensional case, in which each mea- 
surement operator can be a two- or a three-dimensional 
projector, the four-dimensional case with one-, two- 
and three-dimensional projectors, and the three- and 
two-dimensional cases with projectors of any dimen- 
sionality, including degenerate operators, i.e., zero and 
unity. The eight-dimensional calculation covers the 
seven-dimensional case with four- and three-dimensional 
projectors, the six-dimensional one with four-, three- 
and two-dimensional projectors, the five-dimensional one 
with four-, three-, two- and one-dimensional projectors, 
and the four-, three- and two-dimensional ones with any 
projectors. 

The Schmidt decomposition of the maximally violat- 
ing state vector tells which subspaces of the component 
spaces it occupies. The dimensionality of those subspaces 
is the Schmidt number, the number of nonzero Schmidt 
coefficients. If this is less than the number of dimen- 
sions of the component spaces, and if each measurement 
operator projects the corresponding subspace onto itself, 
then the orthogonal subspace may be dropped altogether, 
reducing the dimensionality of the solution. The latter 
condition is fulfilled if the matrix of each measurement 
operator is block-diagonal in the Schmidt basis, such that 
all matrix elements connecting the subspace the state 
vector occupies with the orthogonal subspace are zero. 
Although this was not always true for the results we got 
from the optimization procedure, it turned out that we 
could always get rid of the unwanted nonzero matrix ele- 
ments without changing the optimum value achieved by 



penalizing them in a subsequent optimization. In these 
cases there were equally good solutions with different 
choices of the measurement operators concerned. Thus a 
Schmidt number smaller than the dimensionality of the 
component Hilbert space we considered always led to a 
solution in a smaller space. 

The number of dimensions of a measurement operator 
in the reduced space is given by the dimensionality of the 
overlap of the subspace it projects to with the subspace 
the state vector occupies. We can easily get it as the trace 
of the corresponding block of the matrix of the operator. 
The dimensionality of the overlap of a five-dimensional 
and a three-dimensional subspace of a six-dimensional 
space is either two or three. This is the reason our six- 
dimensional calculation with three-dimensional projec- 
tors may lead to a five-dimensional solution with both 
two- and three-dimensional measurement operators. To 
try all these possibilities for a Bell inequality with five 
measurement settings per party with confining ourselves 
to five-dimensional component spaces would mean 2^" 
separate optimizations. To check all potential four- 
dimensional solutions a single six-dimensional calculation 
covers would require 3^° complete calculations. Although 
the number of parameters to be optimized is larger, it is 
still better to work in higher-dimensional spaces. 

In our previous work [l^ we studied many of the Bell 
inequalities we consider here, but then we confined our- 
selves to four-dimensional component spaces with two- 
dimensional projectors. We solved the problem by max- 
imizing the maximum eigenvalue of the 16 x 16 matrix 
of the operator given by Eq. ^ with a simplex uphill 
method. A difficulty with this approach is that it almost 
always finds a solution with qubits, even if a better so- 
lution exists with higher dimensional component spaces. 
For some inequalities it finds the true maximum only once 
in several thousand attempts. With four-dimensional 
spaces the number of calculations necessary to get re- 
liable results were well affordable. However, in six- or 
eight-dimensional component spaces it is even more diffi- 
cult to find the interesting solution, there are a lot more 
parameters to be optimized, and the size of the matrices 
to be diagonalized many times in each optimization run is 
also much larger. Therefore, we have changed our strat- 
egy and we have used an iterative approach. The proce- 
dure ensures that the Schmidt bases of the state vector 
agree with the bases used in the component spaces, up 
to possible phase factors. 

When we try to find the best m-dimensional solu- 
tion with the n-dimensional (n = 6 or 8) program, we 
first take the rn-dimcnsional maximally entangled state 
l/y^^j^]^ ^ I We apply the simplex uphill method 
[3^ to find the best measurement operators to maxi- 
mize the expectation value of the operator ([4]), that is 
the quantum value of the inequality with the state cho- 
sen. Then we replace the vector with the one of the 
form X]i=i n that gives the largest violation with 

the measurement operators we have got. This vector can 
directly be determined by finding the eigenvector belong- 
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ing to the largest eigenvalue of a matrix. We repeat these 
steps until we reach convergence. As the vector always 
has the special form in our basis, the matrices whose ex- 
pectation value has to be calculated and the ones to be 
diagonalized are of the size of ti x n. We only have to deal 
with small submatrices of the corresponding matrices of 
size X in the full Hilbcrt space. For the sirnplex up- 
hill method we use the implementation of Ref . [33| . The 
procedure does not ensure to find the true optimum, so 
we have to repeat it many times. To save computation 
time we made just a few iterations with many different 
initial values, and chose only a few of them that arrived 
at the best values for the full calculation. In another im- 
plementation we stopped the iterations as soon as any of 
the Schmidt coefficients (which are just \ci\) fell below a 
threshold, and restarted the procedure with new initial 
values. 

To apply optimization we have to parametrize the 
operators. An (ri/2)-dimensional projector in an n- 
dimensional space requires real parameters to 

characterize in the real case, and twice as many in the 
complex case. In choosing the parameters we used a 
similar approach than in Ref. In Ref. [U it has 

been shown that any (n/2)-dimensional projector in an 
n-dimensional space can be transformed into a simple 
form with two unitary transformations, one acting only 
on the first n/2, while the other on the other n/2 basis 
vectors. This simple form is characterized by n/2 angles. 
Then any such projector can be produced by acting on 
such a simple matrix with two appropriate unitary — or 
in the real case orthogonal — transformations. A general 
three-dimensional orthogonal transformation can be de- 
scribed with 3 Euler angles, while the four-dimensional 
one requires 6 Euler angles (see for example Ref. [ssj). 
This way we get the correct number of parameters for the 
projectors in the real case both in the six- (3 -I- 2 • 3 = 9) 
and in the eight-dimensional (4-1-2-6 = 16) spaces. 
In a six-dimensional space each three-dimensional uni- 
tary transformation requires further 4 phase factors to be 
characterized, but one overall factor is irrelevant. The 9 
extra parameters we get this way give just the necessary 
number of parameters for the projectors in the complex 
case. 



IV. UPPER BOUNDS 

To determine upper limits for the maximum quantum 
violation of the Bell inequalities we follow the procedure 
given by Navascues et al. [l^, [13]. The method, in its 
simplest form, specified to the case of a two-party two- 
outcome Bell inequality, is as follows. 

Let us introduce the operators a; = {2Ai — Ia) 'S) Ib 
and bj = Ia® C^Bj — Ib), where I a and Ib are the unity 
operators in Ahce's and Bob's component spaces, respec- 
tively. The Qi and bj operators have eigenvalues -1-1 and 
— 1 for the two outcomes of the measurement instead of 
1 and 0. The square of each of them is unity, and ev- 



ery Ui commutes with every bj. Let us choose a set of 
vectors containing IV)), ai|V')i 02!^'), ■ • ■ ,amAlV')i 
62 . . . , b„is \tp), and any number of vectors of the form 
aiaiiQi'i ...bjbjibjii\ijj), i.e., lip) multiplied by the product 
of any numbers of a and b operators. Let us consider 
the matrix whose entries are the real parts of the scalar 
products of these vectors. A matrix derived this way is 
positive semidcfinite. Each entry of the matrix is the real 
part of the expectation value of an operator. Taking into 
account that each and bj squares to unity, and that 
commutes with bj, one can see that in some places the 
operators may be equal, or hermitian conjugates of each 
other. In those places the entries of the matrix will be 
equal, independently of the actual choice of \tjj) and the 
measurement operators. It is easy to sec that the quan- 
tum value for a Bell inequality can be written as a linear 
combination of certain entries of this matrix. There are 
limits on the value this linear combination may take, as 
the matrix has to be positive semidcfinite, while satisfy- 
ing constraints requiring the equality of some of its ele- 
ments. The quantum value of the Bell inequality can not 
be larger than the maximum value allowed by the con- 
ditions above. The standard optimization method to de- 
termine this maximum is semidcfinite programming [36| . 
There are several codes available to solve the problem, in 
most of our calculations we used CSDP [s^] . 

The upper limit for the maximum violation of the in- 
equality we get this way depends on how we choose the 
set of vectors we generate the matrix from. Refs. [l^ [l^l 
introduce a systematic hierarchy of choices, where the 
higher levels give ever stricter limits until the exact 
maximum violation is achieved. The first level includes 
only the vectors \ip), ailtp) [i ~ l,...,mA) and bj\tp) 
{j = 1, . . . ,mB)- The second level also includes vec- 
tors aitti'lip), bjbjiltjj) and aibj\ip) = 1,...,toa, 
~ 1, . . . tTUb), that is it includes all vectors we can 
make by acting on lip) with up to two measurement oper- 
ators. We get level three by going up to three operators, 
and so on. The exact value for the maximum violation 
may be reached at a finite level, for correlation type Bell 
inequalities this already happens at the first level [38| . It 
has been shown that the series of upper limits through the 
hierarchy converges to the exact value [l^, Slj . In prac- 
tice, it is also worth considering intermediate levels. The 
computation requirement is often much less than that is 
for the next full level, and if it gives the exact value, we 
need not go any further. For example, for many of the 
Bell inequalities we considered we got the upper limit 
equal to the lower one by supplementing level two by 
just the vectors generated with two a and one b (noted 
as aa'b), or one a and two b operators (abb'). 



V. RESULTS FOR THE MAXIMUM 
VIOLATION 

Tableslllland llllj show our results for the maximum vio- 
lation of the Bell inequalities we could achieve by explic- 
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TABLE II: Maximum violation of Bell inequalities. The identitity Ti of the component Hilbert spaces, the dimensionalities Da^ 
and Dsj of the measurement operators of the participants, and the level L when the upper limit reaches the value of violation 
is also shown. Notation la, 2a and 2b mean level 1 plus ab, level 2 plus aa'h (or abb') and level 2 plus aa'b+abb' , respectively. 
An alternative notation for J3322, ^3422 a-nd /4422 is ^43, A4 and A7, respectively. 





Violation 


n 






T 

L 




Violation 
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Db, 


T 

Li 
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n 
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0.2500000 




111 
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no 


Aid 


0.4571068 




0111 


11111 


2a 


A53 


0.6386102 


R^ 


11111 
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2a 


7-1 
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Oil 


nil 


la 


An 


r\ ^ r' A A^^ 

0.3754473 


R^ 


nil 


01111 


2 
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R^ 


11111 
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2a 


7-2 
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R^ 


111 


0111 


2a 


Ais 


0.3843551 


R^ 


nil 


01111 


la 


A 

A55 
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R^ 


01111 
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la 


7-3 

-'3422 


0.4364917 


R^ 
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nil 


la 


AiQ 


0.6226300 




nil 
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2a 


A 
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R^ 


11201 


11111 


2a 


A 

A5 
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R 
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nil 


2a 


A 
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nil 


11111 


2a 


A 
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11111 


11111 


2a 


Ae 


0.3003638 


m)4 

R 
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3 
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0.3258362 


R^ 


2242 


22323 


no 
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0.6488905 


R^ 
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2a 


A C 


0.5412415 


R^ 


1111 


nil 


1 


A 

A22 


0.6234571 




nil 


11111 


2a 


A 

As9 


0.4488256 


R^ 


11111 


01111 


2a 


AS2 


0.8784928 


R 


1111 


nil 


1 


A 

A23 


0.5460735 


/n2 

<L 


11111 


11111 


2 


A 

Aeo 


0.3940032 


R^ 


2224:2 


22243 


2b 


AIIi 


0.6055543 


R^ 


1111 


nil 


la 


A24 


0.6047986 




11111 


11111 


la 


Aai 


0.4019248 




11111 


01111 


2a 
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AII2 


0.5000000 


R^ 


1111 


nil 


2 


A 

A25 


0.6033789 


R^ 


11111 


11111 


la 


A 

A(,2 


0.4043897 


R^ 


22S.i2 


13233 


no 


7-1 

-'4422 


0.2878683 


R'' 


1111 


nil 


2a 


A 

A26 


0.5275550 




11111 


11111 


2 


A 

Afi3 


0.4894164 


R* 


11222 


11222 


2a 


7-2 

-(4422 


0.6213712 


R^ 


1111 


nil 


2 


A27 


0.6483073 


R^ 


11111 


11111 


2 


A(,4 


0.3900890 


R'* 


12122 


11222 


no 


7-3 

^44,22 


0.4142136 


R^ 


2111 


nil 


la 


A28 


0.6403143 


R^ 


11111 
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la 


^65 


0.3688702 


R^ 
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no 


7-4 
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R^ 
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nil 


2a 


A 
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2a 
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0.4364917 


R 
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nil 


la 


A 
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0.5698209 


R 
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la 


A 
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R^ 
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nil 


la 
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R^ 
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la 
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R^ 
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23332 


no 


7-7 
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0.4548373 




1111 


nil 


2 


A32 


0.4135530 




11111 
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2a 




0.6096103 


R^ 


01111 


11111 


2a 


7-g 

-'4422 


0.4877681 


R'' 


2121 


2112 


2b 


A33 


0.6226313 


m2 

<L 


11111 


11111 


2a 


A 

^70 


0.6052228 


m)2 

R 


01111 


11111 


2 


7-9 

-(4422 


0.4616842 


R^ 


1111 


nil 


2a 


A34 


0.5350117 
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11111 


2a 
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f\ A A f\f\'\ 

0.4490163 


R^ 


01111 
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2a 


rlO 
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0.6139456 


R^ 


1111 


nil 


2a 


A35 


0.6249079 


R^ 


11111 


11111 


2a 


^72 


0.6962822 


R^ 


11111 


11111 


2a 


7-11 

-'4422 


0.6383543 


R^ 


1111 


nil 


la 


A 

A36 


0.4388685 




11111 


11111 


2a 


A 

A73 


0.8831381 


R^ 


11111 


11011 


2a 


rl2 

-'4422 


0.6188142 


R^ 


1111 


nil 


2a 


A37 


0.4868868 
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11111 


2a 
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0.6890694 


R^ 


11111 
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2a 


ri3 

-(4422 


0.4348553 


R^ 
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nil 


2a 


A 

A3S 


0.4699126 
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2a 
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2a 
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0.4794103 


R^ 
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2a 


A 
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0.6172035 
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2 


A 
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11111 


2a 


7-15 
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0.4348553 
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R 
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nil 


2a 


A 
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0.6078638 


R 
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la 
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R 
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2a 
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0.4142136 


R^ 
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la 
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0.4785634 
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2a 
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R^ 
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2 
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R^ 
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nil 


2a 


A 
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0.6198655 


R^ 
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la 
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0.6243153 
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2a 
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0.6429670 


R'^ 
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2a 


A 
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0.6107654 


R 
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la 


Aso 


0.3769863 


R 
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no 
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0.4971707 


R^ 
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2b 


A44 


0.5364942 


R^ 
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2a 
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2a 
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11111 
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0.4608544 




11111 
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R* 
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Aio 
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R* 
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0.8004425 
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itly determining the state vector and the measurement 
operators. The size of the component Hilbert space and 
the dimcnsionahties of the operators are also shown. In 
our earher work we have determined the maximum 
violation for all cases shown in Table |TT] with compo- 
nent Hilbert spaces of up to four dimensions taking two- 
dimensional projectors as measurement operators. For 



the inequalities we follow the notations we used there, 
which are the same as the ones in Rcfs. [ll|, By 
extending our calculations to six-dimensional complex 
and real, and eight-dimensional real component spaces 
we could achieve a larger maximum violation for 19 
of the 112 inequalities. We found 9 cases where four- 
dimensional real spaces were necessary for maximum vi- 
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TABLE III; Maximum violation of Bell inequalities introduced in the present paper The identitity Ti. of the component Hilbert 
spaces, the dimensionalities Da^ and Db^ of the measurement operators of the participants, and the level L when the upper 
limit reaches the value of violation is also shown. Notation la, 2a and 2b mean level 1 plus ab, level 2 plus aa'b (or abb') and 
level 2 plus aa'b+abb' , respectively. 
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0.9716836 




nil 


nil 


2a 


789 
■'4422 


1.0035502 




nil 


nil 


la 


t4 

^4422 


0.6861985 


^2 
C 


1111 


nil 


2a 


747 
>'4422 


0.7644749 


R^ 


nil 


nil 


2a 


790 
■'4422 


0.8398157 




nil 


nil 


la 


t5 

J4422 


0.6520359 




1011 


nil 


2a 


748 
>'4422 


0.7510516 


R® 


3334 


3433 


no 


791 
■'4422 


1.2992769 




nil 


nil 


la 


J4422 


0.4461508 




1011 


nil 


2a 


749 
■^4422 


0.8156923 


R 


3222 


2232 


2a 


792 

■'4422 


1.0648162 


R^ 


nil 


nil 


2a 


t7 

J4422 


0.6057263 


R^ 


nil 


1211 


2a 


750 

>'4422 


0.8555932 


R^ 


nil 


nil 


la 


793 
■'4422 


0.9627261 


R^ 


nil 


nil 


la 


jS 

J4422 


0.4772413 


R^ 


1211 


1121 


2a 


751 

■'4422 


0.6750341 


R^ 


1121 


nil 


2a 


794 

■'4422 


0.5944550 


R^ 


0111 


nil 


2a 


J4422 


0.5000000 


R^ 


nil 


1011 


2 


JaA22 


1.0999033 


R^ 


nil 


nil 


la 


795 
■^4422 


1.0014114 


R^ 


0111 


nil 


2a 


tIO 
J4422 


0.6938425 


R^ 


2222 


2121 


2b 


J4422 


0.8093253 


R 


0111 


1101 


2a 


796 

■'4422 


0.9416515 


C 


nil 


nil 


2a 


7II 
J4422 


0.4684758 


r" 


2222 


2121 


2b 


754 
■'4422 


0.5194497 


R^ 


2121 


1221 


2a 


797 

■'4422 


1.1583626 


R^ 


nil 


nil 


no 


t12 
"^4422 


0.7261971 


r" 


2322 


2322 


2b 


755 
>'4422 


0.8283881 


R^ 


2111 


1121 


2a 


798 

■^4422 


1.1496911 


R^ 


nil 


nil 


la 


t13 
J4422 


0.6926870 


R^ 


1121 


nil 


2b 


756 
■'4422 


0.8192037 


R^ 


2111 


1121 


2a 


799 

■^4422 


0.8990726 


R^ 


0111 


nil 


2a 


t14 
J4422 


0.6317469 


R^ 


nil 


1211 


2a 


757 
■'4422 


0.8609885 


R^ 


0111 


1101 


2a 


7IOO 

■'4422 


1.2675376 


R^ 


nil 


nil 


2 


t15 
•^4422 


0.7808500 


R^ 


nil 


1211 


2 


758 
■'4422 


0.8814481 


R^ 


nil 


nil 


la 


7IOI 
■'4422 


1.0296105 


R^ 


nil 


nil 


2a 


t16 
J4422 


0.5674133 


R^ 


nil 


nil 


2a 


759 

■'4422 


0.6379629 


R^ 


0111 


1101 


2 


7IO2 
■^4422 


0.6650865 


R^ 


nil 


nil 


2 


t17 
J4422 


0.6379629 


R^ 


2111 


1121 


2a 


76O 
■'4422 


0.5922714 


R^ 


nil 


nil 


la 


7IO3 
■'4422 


1.0520136 


R^ 


nil 


nil 


2 


t18 
J4422 


1.0130080 


R^ 


1101 


1101 


la 


761 
■'4422 


0.8175538 


R^ 


nil 


nil 


2a 


7IO4 
■^4422 


1.5875885 


R^ 


nil 


nil 


la 


t19 
J4422 


0.6742346 


R^ 


nil 


nil 


la 


762 
■'4422 


0.7500754 


R® 


3333 


3333 


no 


7IO5 
■^4422 


1.0742228 


R^ 


nil 


nil 


la 


t20 
J4422 


0.5932470 


mi2 
R 


1110 


1121 


2a 


763 
■^4422 


0.6078543 


R"' 


2222 


2212 


2a 


7IO6 
■^4422 


0.8382498 


R 


nil 


nil 


la 


t21 
J4422 


0.7017468 


R'^ 


2221 


2212 


2a 


764 

■'4422 


0.9167082 


^2 
C 


nil 


nil 


2a 


7IO7 

■^4422 


0.9339703 


R^ 


nil 


nil 


la 


t22 
J4422 


0.8156179 


R^ 


nil 


nil 


la 


765 
■'4422 


1.1110314 


R^ 


nil 


nil 


2 


7IO8 
■'4422 


0.9676411 


R^ 


nil 


nil 


la 


t23 
J4422 


0.6090289 




nil 


nil 


2a 


766 
■'4422 


0.6186278 


R^ 


nil 


nil 


2a 


7IO9 
■^4422 


1.7260512 


R^ 


nil 


nil 


la 


t24 
J4422 


0.6942505 


R^ 


nil 


nil 


2a 


767 

■'4422 


0.7212948 


R^ 


nil 


1110 


2a 


7IIO 

■^4422 


0.9457465 


R^ 


nil 


nil 


la 


t25 
J4422 


0.5147829 


R^ 


nil 


1211 


2a 


768 
■'4422 


1.0178506 


R^ 


nil 


nil 


la 


7III 

■'4422 


0.7575890 


R^ 


nil 


nil 


la 


t26 
J4422 


0.6402494 


R^ 


nil 


nil 


2a 


769 

■'4422 


0.5081419 


R^ 


nil 


1211 


2a 


7II2 
■'4422 


0.6247585 


R^ 


nil 


nil 


2a 


t27 
J4422 


0.9642857 


^2 
C 


nil 


nil 


2 


770 
■^4422 


0.9538824 


R^ 


nil 


nil 


la 


7I13 
■^4422 


0.8483764 


R^ 


1112 


nil 


3 


t28 
J4422 


0.7500000 


R^ 


nil 


nil 


2a 


771 

■^4422 


0.5822872 


R^ 


nil 


nil 


2a 


7I14 

■'4422 


1.0653155 




nil 


nil 


2a 


7-29 
"^4422 


1.0245760 


^2 
C 


nil 


nil 


2a 


j72 

■'4422 


0.7878251 


C 


nil 


nil 


2a 


7II5 
■^4422 


1.0098044 


R^ 


nil 


1101 


2a 


t30 
J4422 


0.4361950 


R^ 


3231 


2323 


no 


773 

■'4422 


1.1204592 


^2 

c 


nil 


nil 


no 


7II6 

■'4422 


0.7651941 




nil 


nil 


2a 


t31 
J4422 


0.4772413 


R^ 


1211 


1121 


2a 


774 

■'4422 


0.4348553 


R^ 


2111 


nil 


2a 


7II7 

■'4422 


0.9721185 




nil 


nil 


2a 


t32 
J4422 


0.5901592 


R^ 


nil 


nil 


2a 


775 
■'4422 


0.6270262 


R^ 


nil 


nil 


la 


7II8 
■'4422 


0.5283248 




nil 


nil 


2a 


t33 
J4422 


0.6150642 


R^ 


nil 


nil 


2a 


776 

■'4422 


0.8291077 


R^ 


nil 


nil 


la 


7II9 

■'4422 


0.8986088 


R^ 


nil 


1112 


2a 


t34 
J4422 


0.4492657 


R"* 


2222 


2312 


no 


777 
■^4422 


0.8406242 


R2 


nil 


nil 


2 


7I2O 
■'4422 


0.9333181 


R2 


1121 


nil 


2a 


t35 
J4422 


0.6932503 


R2 


nil 


nil 


2a 


778 
■^4422 


0.8937709 


R2 


nil 


1110 


2a 


7I2I 

■^4422 


0.5970802 


R2 


nil 


nil 


2a 


t36 
J4422 


0.6706462 


R2 


nil 


nil 


2a 


t79 

■'4422 


1.0156152 


R2 


nil 


nil 


2a 


7I22 
■'4422 


0.7931031 


R2 


nil 


1101 


2 


t37 
J4422 


0.4864826 


R^ 


1121 


1112 


2a 


78O 

■^4422 


0.9051257 


R2 


nil 


nil 


2 


7I23 
■'4422 


1.0640930 




nil 


nil 


2a 


t38 
J4422 


0.9295153 


R2 


nil 


nil 


2a 


78I 

■'4422 


0.5995871 


R= 


nil 


nil 


2 


7I24 

■'4422 


0.9409648 


9? 


nil 


nil 


2 


739 
J4422 


0.8264535 


R^ 


nil 


nil 


2a 


782 
■'4422 


0.7308078 




nil 


nil 


2a 


7I25 
■'4422 


1.0000000 


R2 


nil 


nil 


2a 


740 
J4422 


0.9208907 


R-* 


2212 


2221 


3 


783 
■'4422 


0.6691786 


R2 


nil 


nil 


2a 


7I26 
■^4422 


0.8023577 




nil 


nil 


2a 


t41 
J4422 


0.7596044 


R2 


nil 


nil 


2 


784 
■'4422 


1.0669560 


R2 


nil 


nil 


2a 


7I27 

■'4422 


0.6718972 


R^ 


1222 


2121 


2a 


t42 
J4422 


0.6722257 


R2 


nil 


1121 


2a 


785 
■'4422 


0.9762793 


R2 


nil 


nil 


2 


7I28 
■'4422 


1.0096066 


R2 


nil 


nil 


la 


743 
J4422 


0.4685939 


R^ 


2221 


2212 


2a 


786 
■^4422 


0.7558762 


R^ 


1112 


nil 


2b 


7I29 
■'4422 


1.0522569 


R2 


nil 


nil 


la 
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TABLE IV: Difference between the upper limit calculated at different levels and the lower limit derived by explicitly constructing 
the measurement operators. 



Case 


L2 


L2 + aa'h 


L2 + 85% of L2 + aa'b+ 


L2 + aa'b+ 


L3 


L3 + aa'a"b+ 








aa'b + abb' 


+abb' 


abb' + aa'a" 




aa'bb' + abb'b" 


^3322 


0.0009397 


0.0008863 




0.0008760 


0.0008756 


0.0008756 


0.0008754 


t20 
J 4422 


0.0393814 


0.0111147 




0.0040799 




0.0000784 




t30 
J4422 


0.0401921 


0.0109109 




0.0042758 




0.0010680 




734 
"^4422 


0.0240714 


0.0064382 




0.0023265 




0.0001503 




t48 
J4422 


0.0294613 


0.0014661 




0.0008068 




0.0006361 




t62 
J4422 


0.0022695 


0.0000194 




0.0000178 




0.0000167 




t73 
J4422 


0.0329957 


0.0030318 




0.0023042 




0.0020987 




t97 
J4422 


0.0392054 


0.0018759 




0.0005867 




0.0001425 




Ai4 


0.0274655 


0.0057079 




0.0031815 


0.0025615 






A21 


0.0048537 


0.0007595 




0.0005165 


0.0003840 






Ai7 


0.0202626 


0.0050055 


0.0026688 










Ag2 


0.0133913 


0.0031907 


0.0025556 










AliA 


0.0164779 


0.0032730 


0.0016572 










^65 


0.0175447 


0.0014105 


0.0000806 










AdJ 


0.0135307 


0.0019890 


0.0006118 










Ads. 


0.0073408 


0.0039257 


0.0038884 












0.0161885 


0.0026707 


0.0006052 










^82 


0.0082846 


0.0065568 


0.0064928 










^84 


0.0271831 


0.0040203 


0.0021524 










^89 


0.0104965 


0.0043258 


0.0033485 











olation, with not all projectors two-dimensional. For an- 
other 9 inequalities we had to go up to five-dimensional 
real component spaces. One of them is ^10, a 4522 in- 
equality, which has less than five measurement settings 
for one of the parties. For one 5522 we found the max- 
imum violation with six-dimensional component spaces, 
more than the number of measurement settings for either 
party. Similarly, for J4422 (see Table HIT)), we needed five- 
dimensional, and for J4422 ^^'^ ^^44221 4422 cases, six- 
dimensional real component spaces. However, for these 
cases we can not be sure that the states and the operators 
we constructed does give the true quantum bound. The 
values for the violation arc still somewhat below the cor- 
responding upper limits we have got. Some possibilities, 
even with lower dimensional component spaces, like four- 
dimensional complex spaces with allowing degenerate 
measurement operators are not covered. Moreover, we 
may have even missed some solutions which are in prin- 
ciple accessible by our programs. The six-dimensional 
program with complex spaces often struggles to find so- 
lutions it finds relatively easily in real spaces. This is not 
very surprising, as in the complex case there are twice as 
many parameters to optimize. We note that whenever we 
needed more than two-dimensional component spaces for 
maximum violation for the inequalities considered in the 
present paper, real spaces were always sufficient. This 
is not true for all Bell inequalities, for example for the 
4822 correlation type inequality introduced in Ref. [i^ 



the lowest dimensional component spaces to achieve the 
maximum violation are four-dimensional complex ones. 
Some complex solutions may have been missed by our 
present calculations, but not many of them, as for the 
great majority of the inequalities the solution wc man- 
aged to construct are actually as good as they can be, 
according to the upper limit we found. 

For the cases the upper limit agreed with the lower one 
derived from the explicit construction, we also show in 
Tables and IIIII the (usually partial) level when it hap- 
pened. For those inequalities the maximum violations 
shown arc the exact values. The level wc could afford 
to accomplish with the code and with the computers we 
used depended on the number of measurement settings. 
For the 5522 inequalities we could do the calculation with 
either all aa'b or abb' type vectors added to level two, but 
with not both of them. Therefore, the largest calculations 
we made included only the 85% of them, randomly cho- 
sen. For the 4522 cases the maximum we could do was 
level two plus aa'b + abb' + aa'a" (i.e., level 3 without 
bb'b"). For 4422 wc could go up to level 3, while for 3322 
to level three plus aa' bb" + aa' a" b + abb' b" . Wc had to go 
to the maximum we could afford only for a small fraction 
of the inequalities. For over 85% of the cases level two 
plus abb' (and for the same cases level two plus aa'b) has 
already given the exact value. 

There were 20 inequalities for which the upper limit 
remained above the best solution we found (marked by 
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no in the Tables) at all levels we could afford. The gap 
between the upper and lower limits for those inequalities, 
depending on the level is shown in Table Hvl In almost 
all of these cases extending the calculation for the higher 
limit further decreases this gap, often quite significantly. 

For /|^22 J "^4422 1 "^4422 1 "^4422 J ^65 , ^67 and Ago it sccms 
to be probable that the lower limits we got are exact. 
When the decrease of the gap is relatively very small, 
like in the cases of J4422, ^68 or ^§2, it is more hkely 
that we still have not got the true solution. 

However, we find the case of 13322 truly puzzling. This 
is one of the simplest tight Bell inequalities, with just 
three measurement settings for each party. In our formu- 
lation used in ^ for the Bell inequalities it reads 

-2pAi ~PA2 -PBi +PA1B1 +PA1B2 +PA1B3 

+ PA2B1 +PA2B2 -PA2B3 +PA3B1 -PA3B2 < 0. (5) 

We could calculate the upper limit at a level significantly 
exceeding level 3 (up to level 3 the upper limit has also 
been calculated by Refs. [10, HH). The gap between the 
lower limit of 0.25 and the upper limit is at least 4 orders 
of magnitude larger than numerical uncertainty, while 
it gets just marginally smaller by the extension of the 
calculation for the upper limit. It does not behave as 
if it would decrease much further. The lower limit has 
already been achieved with real two-dimensional compo- 
nent spaces. All our attempts to find a larger violation 
with higher dimensional spaces failed, although the num- 
ber of parameters is much smaller than for the other in- 
equalities we considered, and we have also made a lot 
more attempts than for any other case. It does not seem 
very likely for us that with spaces and operators our cal- 
culations cover such a solution exists. At the same time, 
it would also be very surprising, if one needed even higher 
dimensional Hilbert spaces to violate maximally this very 
simple inequality. 

Our present results confirm what we have found in 
Ref . [H I , namely in most cases the state giving the max- 
imum violation is not the maximally entangled one (in- 
dication for this fact was also given independently by 
Ref. [l3|)- The results here are stronger, because the 
comparison with the upper bound proves that for most 
inequalities the states and operators we have got are actu- 
ally the ones giving the absolute quantum bound. When 
the component Hilbert spaces required were larger than 
two-dimensional, the state we found was never the max- 
imally entangled one. For bipartite inequalities with two 
measurement settings, but more than two outcomes (39| . 
a similar result was found by Ref. [ll] . 



VI. DETECTION EFFICIENCIES 

A loophole- free experiment in order to test the nonlocal 
nature of quantum mechanics is still missing 01. In par- 
ticular, none of the experiments performed to date could 
close simultaneously the locality loophole (the measure- 



ment results at Alice's and Bob's side should be space- 
like separated) and the detection loophole. In order to 
avoid this latter loophole the particles must be detected 
with a high enough probability, otherwise a locally causal 
model can reproduce the measured correlations. Another 
motivation beside the fundamental ones comes from the 
security issue of some quantum communication protocols 
which is based on the loophole-free violation of Bell in- 
equalities [4^ . 

There are different proposals to close the detection 
loophole. More than two outcomes [H, I3, more than 
"MJil S H m, and partially entangled 
46l |47I has been considered as well. The case 



two settings 
states [IT' 

of asymmetric Bell experiments, where the two particles 
are detected with different probabilities (in systems such 
as entaiigled atom-photon pairs) has also been addressed 
[13, [13, [43 ■ Another promising approach is the applica- 
tion of homodyne detectors in Bell tests [H, [4^, [s^ • 

Here we calculate both the symmetric (Alice and Bob 
have the same efficiency tja ^ Vb ^ v) and asymmetric 
(Alice's detector is perfect) threshold detection efficien- 
cies of our set of Bell inequalities for a pair of maximally 
entangled qubits allowing degenerate measurements as 
well. The general approach, consisting of two different 
detection efficiencies {rjA and tje) was treated in Ref. [II]- 
According to their Eq. (4) the quantum value of a Bell 
inequality with detection efficiencies rjA and rjB is given 
by 

Iru-VB = iUVbQ + ??a(1 - ■riB)MA 

+ (1 - iia)vbMb + (I - iia){1 - Vb)X, (6) 

where Q is the quantum value associated to the Bell in- 
equality with perfect detectors, Ma,b and X are the val- 
ues when one or both detectors do not fire. Next we limit 
ourselves to maximally entangled states and in the case 
of no detection one of the two outcomes {1,0} are taken 
as an output. If we consider degenerate measurements as 
well, the no-detection outcomes must be set to the same 
value as the output of the corresponding degenerate mea- 
surement. On the other hand, by fixing these values the 
original Bell inequality / reduces to another one /' with 
a smaller number of settings. The local bound L for this 
inequality can be smaller or equal to the local bound cor- 
responding to the Bell inequality /. At the same time it 
can be observed that the values Ma,b and X remain the 
same for /'. Thus, if for a given Bell inequality / and ef- 
ficiencies r]A,r]B, IriA,riB > ^ holds, it certainly holds true 
for the reduced Bell inequality /' as well. This implies 
that the threshold efficiencies for /' can only be lowered 
with respect to the original one / (in our actual calcu- 
lations we found that L never decreased, hence leaving 
unchanged the detection efficiencies). Conversely, this 
result means that if one supplements a Bell inequality 
with coefficients pertaining to degenerate measurements 
(i.e., deterministic ones, which need not be performed at 
all), it cannot lower the detection efficiency. Note, that 
this result, which is intuitively clear, also holds for the 
case of partially entangled states. We made an optimiza- 
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tion over all possible measurement strategies (i.e., for any 
combination of degenerate and non-degenerate measure- 
ments) and over all possible no-detection outcomes, both 
in the symmetric and asymmetric situations. The follow- 
ing results have been obtained: 

From the literature [l^] for the maximally entangled, 
asymmetric case the best known detection efficiency is 
rjB = 2/3, corresponding to 13322- We found that in- 
equality /I422 and AII2 yield the same value with non- 
degenerate measurements. With involving degenerate 
measurements as well, this threshold value was obtained 
for a further 97 inequalities. Out of them 69 can be 
traced back to I3322, while 12 is reducible to non-tight 
3322 inequalities, nonequivalent with 13322. However, 
T]B = 2/3 is not the best one can achieve, for three 
five-setting Bell inequalities from the list of Avis et al. 
[ll| . namely, A34, A44, and Ac,o do better. The respec- 
tive values of rjs = 0.6607, 0.6520, 0.6587 are reached 
with purely non-degenerate measurements. There exists 
a local model [5l| , which reproduces the correlations aris- 
ing from non-degenerate projective measurements which 
can be expressed using real (complex) numbers on the 
maximally entangled state under the assumption that 
Alice has perfect detector and Bob has detection effi- 
ciency r]B ~ 2/3 {i]B = 1/2). Thus the three lower than 
2/3 values given above need to belong to projective mea- 
surements with settings requiring complex numbers to be 
described. This fact, as we have checked, is true. 

For the maximally entangled symmetric case, the 
threshold efficiency for the CHSH inequality is known 
to be 0.8284, which result was slightly improved recently 
by Brunner and Gisin l^. They showed that the A5 



inequality from the Avis et al. list allows the slightly 
smaller threshold of 0.8214. By computing the best 
threshold values for the set of 241 Bell inequalities we 
did not find any better case. Allowing degenerate mea- 
surements we could get the same value for 19 five-setting 
Bell inequalities, but they could all be traced back to A5 . 



VII. SUMMARY 

Let us summarize the main results achieved in this 
work. 

We have presented two heuristic methods to get tight 
two-party two-outcome Bell inequalities. Using these 
methods we have extended the list of known 4422 type 
inequalities by 129 members. The list is very probably 
still not full. 

For 241 inequalities we used numerical optimization 
to determine measurement operators and states giving 
their maximum quantum violation achievable with six- 
dimensional complex component Hilbert spaces, allow- 
ing three-dimensional projectors as measurement opera- 
tors, and with eight-dimensional real component Hilbert 
spaces, allowing four-dimensional projectors as measure- 



ment operators. An n-dimensional calculation with 
(ri/2)-dimcnsional projectors covers all m-diniensional 
(to < n) cases with projectors having any dimension- 
ality between max(0, to — n/2) and min(m, n/2). The 
241 examples we have considered include all tight bipar- 
tite Bell inequalities with up to five two-outcome mea- 
surement settings known to us (excluding CHSH, about 
which we could tell nothing new). These results represent 
lower limits for the quantum bounds for the maximum 
violation. At the same time, the method proposed in 
Refs. [1^ [lO] makes it possible to determine a series of 
upper bounds by carrying out the calculation through 
a hierarchy of levels. The series converges to the exact 
value, often reaching it at a finite level. For the majority 
of cases the upper bound has become equal to the lower 
one at a level which is quite easy to perform for inequal- 
ities with no more than five settings per party. There 
remained only 20 cases for which we still can not tell for 
sure the maximum value of their quantum violation. The 
lower limit is probably the exact value for quite a few of 
them, but unfortunately we could not do the calculation 
for the upper limit at a high enough level to verify it. It is 
surprising that for the smallest case we considered, with 
only three measurement settings per party, there is still 
a significant gap between the upper and lower bounds, 
while the former seems to come only marginally lower by 
increasing the level, and for the latter we were unable to 
find any better solution than one can get with real qubits, 
despite the high dimensionality of the component spaces 
we considered and the many attempts we have done. We 
have only found a few inequalities whose maximum vio- 
lation was achieved with the maximally entangled state, 
in all those cases a pair of qubits were sufficient. 

We have also calculated the minimum detection effi- 
ciency to verify quantum violation for each inequality 
with a maximally entangled pair of particles, both for 
symmetric and asymmetric arrangements. By allowing 
degenerate measurements (ones with a certain outcome) , 
one could lower these threshold efficiencies for most in- 
equalities, but this only means that the (not necessarily 
tight) inequality with less measurement settings one can 
get by skipping the degenerate ones is at least as good at 
verifying quantum violation. As far as the actual num- 
bers are concerned, for the symmetric case none of the 
inequalities does better than the best one already known, 
while for the asymmetric case we have found three in- 
equalities with five measurement settings per party, for 
which slightly less efficient detectors are sufficient than 
for the best one known so far. 
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